V-primer: software for the efficient design of genome-wide InDel and SNP markers from multi-sample variant call format (VCF) genotyping data.

DNA markers are indispensable tools in genetics and genomics research as well as in crop breeding, particularly for marker-assisted selection. Recent advances in next-generation sequencing technology have made it easier to obtain genome sequences for various crop species, enabling the large-scale identification of DNA polymorphisms among varieties, which in turn has made DNA marker design more accessible. However, existing primer design software is not suitable for designing many types of genome-wide DNA markers from next-generation sequencing data. Here, we describe the development of V-primer, high-throughput software for designing insertion/deletion, cleaved amplified polymorphic sequence, and single-nucleotide polymorphism (SNP) markers. We validated the applicability of these markers in different crops. In addition, we performed multiplex PCR targeted amplicon sequencing using SNP markers designed with V-primer. Our results demonstrate that V-primer facilitates the efficient and accurate design of primers and is thus a useful tool for genetics, genomics, and crop breeding. V-primer is freely available at https://github.com/ncod3/vprimer.

Disentangling the complex gene interaction networks between rice and the blast fungus identifies a new pathogen effector.

Studies focused solely on single organisms can fail to identify the networks underlying host-pathogen gene-for-gene interactions. Here, we integrate genetic analyses of rice (Oryza sativa, host) and rice blast fungus (Magnaporthe oryzae, pathogen) and uncover a new pathogen recognition specificity of the rice nucleotide-binding domain and leucine-rich repeat protein (NLR) immune receptor Pik, which mediates resistance to M. oryzae expressing the avirulence effector gene AVR-Pik. Rice Piks-1, encoded by an allele of Pik-1, recognizes a previously unidentified effector encoded by the M. oryzae avirulence gene AVR-Mgk1, which is found on a mini-chromosome. AVR-Mgk1 has no sequence similarity to known AVR-Pik effectors and is prone to deletion from the mini-chromosome mediated by repeated Inago2 retrotransposon sequences. AVR-Mgk1 is detected by Piks-1 and by other Pik-1 alleles known to recognize AVR-Pik effectors; recognition is mediated by AVR-Mgk1 binding to the integrated heavy metal-associated (HMA) domain of Piks-1 and other Pik-1 alleles. Our findings highlight how complex gene-for-gene interaction networks can be disentangled by applying forward genetics approaches simultaneously to the host and pathogen. We demonstrate dynamic coevolution between an NLR integrated domain and multiple families of effector proteins.

Genome Analysis Revives a Forgotten Hybrid Crop Edo-dokoro in the Genus Dioscorea.

A rhizomatous Dioscorea crop 'Edo-dokoro' was described in old records of Japan, but its botanical identity has not been characterized. We found that Edo-dokoro is still produced by four farmers in Tohoku-machi of the Aomori prefecture, Japan. The rhizomes of Edo-dokoro are a delicacy to the local people and are sold in the markets. Morphological characters of Edo-dokoro suggest its hybrid origin between the two species, Dioscorea tokoro and Dioscorea tenuipes. Genome analysis revealed that Edo-dokoro likely originated by hybridization of a male D. tokoro to a female D. tenuipes, followed by a backcross with a male plant of D. tokoro. Edo-dokoro is a typical minor crop possibly maintained for more than 300 years but now almost forgotten by the public. We hypothesize that there are many such uncharacterized genetic heritages passed over generations by small-scale farmers that await serious scientific investigation for future use and improvement by using modern genomics information.

A genetically linked pair of NLR immune receptors shows contrasting patterns of evolution.

Throughout their evolution, plant nucleotide-binding leucine-rich-repeat receptors (NLRs) have acquired widely divergent unconventional integrated domains that enhance their ability to detect pathogen effectors. However, the functional dynamics that drive the evolution of NLRs with integrated domains (NLR-IDs) remain poorly understood. Here, we reconstructed the evolutionary history of an NLR locus prone to unconventional domain integration and experimentally tested hypotheses about the evolution of NLR-IDs. We show that the rice (Oryza sativa) NLR Pias recognizes the effector AVR-Pias of the blast fungal pathogen Magnaporthe oryzae. Pias consists of a functionally specialized NLR pair, the helper Pias-1 and the sensor Pias-2, that is allelic to the previously characterized Pia pair of NLRs: the helper RGA4 and the sensor RGA5. Remarkably, Pias-2 carries a C-terminal DUF761 domain at a similar position to the heavy metal-associated (HMA) domain of RGA5. Phylogenomic analysis showed that Pias-2/RGA5 sensor NLRs have undergone recurrent genomic recombination within the genus Oryza, resulting in up to six sequence-divergent domain integrations. Allelic NLRs with divergent functions have been maintained transspecies in different Oryza lineages to detect sequence-divergent pathogen effectors. By contrast, Pias-1 has retained its NLR helper activity throughout evolution and is capable of functioning together with the divergent sensor-NLR RGA5 to respond to AVR-Pia. These results suggest that opposite selective forces have driven the evolution of paired NLRs: highly dynamic domain integration events maintained by balancing selection for sensor NLRs, in sharp contrast to purifying selection and functional conservation of immune signaling for helper NLRs.

Sat-BSA: an NGS-based method using local de novo assembly of long reads for rapid identification of genomic structural variations associated with agronomic traits.

Advances in next generation sequencing (NGS)-based methodologies have accelerated the identifications of simple genetic variants such as point mutations and small insertions/deletions (InDels). Structural variants (SVs) including large InDels and rearrangements provide vital sources of genetic diversity for plant breeding. However, their analysis remains a challenge due to their complex nature. Consequently, novel NGS-based approaches are needed to rapidly and accurately identify SVs. Here, we present an NGS-based bulked-segregant analysis (BSA) technique called Sat-BSA (SVs associated with traits) for identifying SVs controlling traits of interest in crops. Sat-BSA targets allele frequencies at all SNP positions to first identify candidate genomic regions associated with a trait, which is then reconstructed by long reads-based local de novo assembly. Finally, the association between SVs, RNA-seq-based gene expression patterns and trait is evaluated for multiple cultivars to narrow down the candidate genes. We applied Sat-BSA to segregating F2 progeny obtained from crosses between turnip cultivars with different tuber colors and successfully isolated two genes harboring SVs that are responsible for tuber phenotypes. The current study demonstrates the utility of Sat-BSA for the identification of SVs associated with traits of interest in species with large and heterozygous genomes.

Genome analyses reveal the hybrid origin of the staple crop white Guinea yam (Dioscorea rotundata).

White Guinea yam (Dioscorea rotundata) is an important staple tuber crop in West Africa. However, its origin remains unclear. In this study, we resequenced 336 accessions of white Guinea yam and compared them with the sequences of wild Dioscorea species using an improved reference genome sequence of D. rotundata In contrast to a previous study suggesting that D. rotundata originated from a subgroup of Dioscorea praehensilis, our results suggest a hybrid origin of white Guinea yam from crosses between the wild rainforest species D. praehensilis and the savannah-adapted species Dioscorea abyssinica We identified a greater genomic contribution from D. abyssinica in the sex chromosome of Guinea yam and extensive introgression around the SWEETIE gene. Our findings point to a complex domestication scenario for Guinea yam and highlight the importance of wild species as gene donors for improving this crop through molecular breeding.

Next-generation sequencing-based bulked segregant analysis for QTL mapping in the heterozygous species Brassica rapa.

KEY MESSAGE: An improved protocol of QTL-seq, an NGS-based method for bulked segregant analysis we previously developed in rice, allowed successful mapping of QTLs of interest in the highly heterozygous genome of B. rapa, demonstrating the power of this elegant method for genetic analyses in heterozygous species of economic importance. Recent advances in next-generation sequencing (NGS) and the various NGS-based methods developed for rapidly identifying candidate genes of interest have accelerated genetic analysis mainly in the model plants rice and Arabidopsis. Brassica rapa includes several economically important crops such as Chinese cabbage, turnip and various leafy vegetables. The application of NGS-based approaches for the analysis of B. rapa has been limited mainly due to its highly heterozygous genome and poor quality of the reference genome sequence currently available for this species. In this study, we have improved QTL-seq, a method for NGS-based bulked segregant analysis we previously developed in rice, extending its applicability for accelerating the genetic analysis and molecular breeding of B. rapa. Addition of new filters to the original QTL-seq pipeline allowed removal of spurious single-nucleotide polymorphisms caused by alignment/sequencing errors and variability between parents, significantly improving accuracy of the analysis. As proof of principle, we successfully applied the new approach to identify candidate genomic regions controlling flowering and trichome formation using segregating F2 progeny obtained from crosses made between cultivars of B. rapa showing contrasting phenotypes for these traits. We strongly believe that the improved QTL-seq method reported here will extend the applicability of NGS-based genetic analysis not only to B. rapa but also to other plant species of economic importance with heterozygous genomes.

Body balance function of cochlear implant patients with and without sound conditions.

OBJECTIVE: The relation between well-controlled auditory stimulation through cochlear implant (CI) and the body balance has been sparsely investigated. The purpose of this study was to evaluate the body balance function of CI patients with- and without-sound in anechoic sound-shielded room. METHODS: We recorded 8 experienced CI recipients and 8 young normal-hearing volunteers. All subjects were assessed using posturography under 4 conditions: (1) eyes open with-sound, (2) eyes closed with-sound, (3) eyes open without-sound, and (4) eyes closed without-sound. RESULTS: The total path length and the total area were significantly larger in the eyes closed condition than in the eyes open condition. In normal hearing subjects, the average displacement of center of pressure (COP) in the mediolateral direction under with-sound condition was not different from that under without-sound condition. In CI recipients, the COP significantly displaced to the CI side after the deprivation of visual cues in without-sound condition. This shift was eliminated in with-sound condition (significant interaction among sound condition, eye condition, and between-group factor). CONCLUSION: In CI subjects, sound stimulation improves the abnormal displacement of COP in the mediolateral direction. SIGNIFICANCE: A posturographic study under an anechoic condition proved that sound stimulation improves body balance function in CI subjects.

Genome sequencing of the staple food crop white Guinea yam enables the development of a molecular marker for sex determination.

BACKGROUND: Root and tuber crops are a major food source in tropical Africa. Among these crops are several species in the monocotyledonous genus Dioscorea collectively known as yam, a staple tuber crop that contributes enormously to the subsistence and socio-cultural lives of millions of people, principally in West and Central Africa. Yam cultivation is constrained by several factors, and yam can be considered a neglected "orphan" crop that would benefit from crop improvement efforts. However, the lack of genetic and genomic tools has impeded the improvement of this staple crop. RESULTS: To accelerate marker-assisted breeding of yam, we performed genome analysis of white Guinea yam (Dioscorea rotundata) and assembled a 594-Mb genome, 76.4% of which was distributed among 21 linkage groups. In total, we predicted 26,198 genes. Phylogenetic analyses with 2381 conserved genes revealed that Dioscorea is a unique lineage of monocotyledons distinct from the Poales (rice), Arecales (palm), and Zingiberales (banana). The entire Dioscorea genus is characterized by the occurrence of separate male and female plants (dioecy), a feature that has limited efficient yam breeding. To infer the genetics of sex determination, we performed whole-genome resequencing of bulked segregants (quantitative trait locus sequencing [QTL-seq]) in F1 progeny segregating for male and female plants and identified a genomic region associated with female heterogametic (male = ZZ, female = ZW) sex determination. We further delineated the W locus and used it to develop a molecular marker for sex identification of Guinea yam plants at the seedling stage. CONCLUSIONS: Guinea yam belongs to a unique and highly differentiated clade of monocotyledons. The genome analyses and sex-linked marker development performed in this study should greatly accelerate marker-assisted breeding of Guinea yam. In addition, our QTL-seq approach can be utilized in genetic studies of other outcrossing crops and organisms with highly heterozygous genomes. Genomic analysis of orphan crops such as yam promotes efforts to improve food security and the sustainability of tropical agriculture.

A chloroplast-localized protein LESION AND LAMINA BENDING affects defence and growth responses in rice.

Understanding how plants allocate their resources to growth or defence is of long-term importance to the development of new and improved varieties of different crops. Using molecular genetics, plant physiology, hormone analysis and Next-Generation Sequencing (NGS)-based transcript profiling, we have isolated and characterized the rice (Oryza sativa) LESION AND LAMINA BENDING (LLB) gene that encodes a chloroplast-targeted putative leucine carboxyl methyltransferase. Loss of LLB function results in reduced growth and yield, hypersensitive response (HR)-like lesions, accumulation of the antimicrobial compounds momilactones and phytocassanes, and constitutive expression of pathogenesis-related genes. Consistent with these defence-associated responses, llb shows enhanced resistance to rice blast (Magnaporthe oryzae) and bacterial blight (Xanthomonas oryzae pv. oryzae). The lesion and resistance phenotypes are likely to be caused by the over-accumulation of jasmonates (JAs) in the llb mutant including the JA precursor 12-oxo-phytodienoic acid. Additionally, llb shows an increased lamina inclination and enhanced early seedling growth due to elevated brassinosteroid (BR) synthesis and/or signalling. These findings show that LLB functions in the chloroplast to either directly or indirectly repress both JA- and BR-mediated responses, revealing a possible mechanism for controlling how plants allocate resources for defence and growth.